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Abstract 


The Kodiak Islands are located approximately 120 to 250 km from the Alaska- Aleutian Trench 
and are within the southern extent of the 1964 Prince William Sound (M w = 9.2) earthquake rupture 
zone. Here we report new campaign GPS results (1993-2001) from northern Kodiak Island. The rate 
and orientation of the horizontal velocities, relative to a fixed North America, range from 25.3 ±1.4 
mm/yr at N32.9°W ± 2.5° to 8.5 ± 1.0 mm/yr at ^9.7^ ± 6.5°. In addition to the northern Kodiak 
data, we analyzed data from three southern Kodiak Island stations. The inlan d stations from both the 
northern and southern networks indicate a counterclockwise rotation of the velocity vectors. These 
results are consistent with the hypothesis that the difference between the Pacific-North American plate 
motion and the orientation of the down going slab would lead to 4-8 mm/yr of left-lateral slip above the 
unlocked, down-dip portion of the main thrust zone. The northern and southern Kodiak geodetic data 
are consistent with a model that includes the viscoelastic response to (1) a downgoing Pacific plate 
interface that is locked at shallow depths, (2) local coseismic slip in the 1964 earthquake, and (3) 
interseismic creep down dip from the seismogenic zone. Based on the pre-1964 and post- 1964 
earthquake history, as well as the pattern of interseismic earthquakes across the plate boundary zone, we 
hypothesize that in southern Kodiak some strain is released in moderate to large earthquakes between 
the occurrencesof gr &a x ea rt h q u ake s likBthe 4964 -event. - 
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Geologic and Tectonic Setting 

The motion of the Pacific plate relative to the North American plate is at a rate of about 57 
mm/yr at N22°W [DeMets et aL, 1994]. The Kodiak Islands span a region approximately 120 to 250 km 
from the Alaska- Aleutian Trench (taken here to be the eastern extent of the 4000 m or 5200 m 
bathymetric contour). The subducting Pacific plate is thought to have a shallow dip (<10°) from the 
Aleutian trench northwestward to 30-40 km below Kodiak Island based on seismicity [Davies et aL, 
1981; Pulpan and Frohlich, 1985] and the geological structure suggested from a deep crustal transect 
[vonHuene et aL, 1999] (Figure 1). From the northern to the southern part of the island the distance 
from the Aleutian trench to the eastern coast of Kodiak Island narrows by about 10% (Figure 1). This 
change probably reflects a slight steepening of the dip of the down going Pacific plate under the 
southern part of the island. The average trend of the trench off northern Kodiak is approximately N60°E 
which gives a trench perpendicular motion more westerly than the relative plate motion vector by about 
8°. The difference in orientation between the trench normal and the relative plate motion has potential 
implications for the crustal deformation patterns [. McCaffrey et aL, 2000]. 

All of the Kodiak Islands lie within the southern extent of the 1964 Prince William Sound (M w = 
9.2) earthquake rupture and aftershock zone. Prior to the 1964 earthquake, in 1900, a M s = 7.7 



personal communication, 2004]. Nishenko and Jacob [1990] suggest a recurrence interval of 
approximately 60 years between periods of increased activity that last as long as 10 years based on large 
and great earthquakes felt on Kodiak Island during the last 200 years; they considered all of the Kodiak 
Islands as one segment of the Alaska- Aleutian subduction zone. Recently, several large earthquakes 
(M w = 7.0, in 1999; M w = 6.5 in 2000; M w = 7.0 in 2001) have occurred in the southern part of the 
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Kodiak Islands as well [Ratchkovski and Hansen, 2001; Hansen and Ratchkovsfd, 2001]. The Harvard 
and U.S. Geological Survey (USGS) Centroid Moment Tensor solutions (CMTs, Dziewonski and 
Woodhouse , 1983; Ekstrom, 1994; http://neiausgs.gov/neis/sopar) for the events are given in Figure 2. 
The focal mechanisms between Nov. 1977 and Feb. 2004 indicate a combination of reverse faulting 
events on northeast-southwest striking planes and left-lateral strike-slip events at roughly trench-parallel 
orientations (Figure 2). These earthquakes, as well as earthquake focal mechanisms from 1964 - 1979, 
are discussed in relation to the 1964 earthquake in Doser et al., [2002]. They found that the earthquakes 
between 1964 and 1979 were characterized by more normal and normal-oblique mechanisms. The 
earthquakes since 1964, as well as historical earthquakes reported in Gilpin [1995], indicate that the 
southern Kodiak Island region has had moderate to large earthquakes more frequently than northern 
Kodiak both prior to, and following, the 1964 earthquake. 

The Kodiak Islands are part of a large subduction complex that comprises the eastern Aleutian 
forearc; the islands form the subaerial part of a broad topographic ridge that includes a Mesozoic- 
Cenozoic accretionary complex (Figure 1). The Kodiak Seamount (24 Ma) shown in Figure 1 is part of 
the Kodiak - Bowie seamounts, a hot spot chain that transects the Gulf of Alaska; it is the oldest 
seamount still exposed prior to subduction. vonHuene et al. [1999] hypothesized that a thickened and 
thermally hot oceanic crust in this region is a more buoyant segment than adjacent crust and may have 
actedasasec ond seismic asperity in -theT964 earthqu ake (the pr imary asperity amFepicenter was in the 
Prince William Sound area). 

In the next section, we report horizontal velocities and uplift rates estimated from GPS 
measurements made between 1993 and 2001 across the northern region of Kodiak. To put these results 
in a regional context we included in our analysis GPS measurements made by Savage et al, [1999] 
across the southern portion of the Kodiak Islands and measurements from the two permanent stations 
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(KODl and KODK). We used these geodetic results to further address the following questions: (1) how 
do the ongoing crustal deformation rates vary as a function of distance from the trench? (2) Are there 
differences in the deformation rates between northern and southern Kodiak? (3) What was the 
magnitude of 1964 coseismic slip near northern and southern Kodiak? (4) What is the seismic history of 
moderate to large earthquakes in northern and southern Kodiak? Are there differences between these 
two regions? (5) Does the rale of horizontal and vertical deformation vary as a function of time since 
the 1964 earthquake? 


Global Positioning System Measurements and Analysis 

To estimate station velocities we have used observations from eight locations on Kodiak Island 
acquired with varying occupation scenarios between 1993 and 2001 (Table 1), along with observations 
from the two Kodiak permanent GPS stations. We surveyed five primary stations in the northern part of 
the island (PAS A, MILB, CLAM, KODL and SKIO) three or four times between 1993 and 2001, 
occupying most of these for 6-8 hours on 1-4 days but KODl continuously for the length of each survey. 
In 1995, 1997, and 1999 these observations were made as part of an education outreach program with 
Kodiak Island High School [Stockman et aL, 1997; Sauber et aL, 1998; E. Linscheid 
http://137.229.234.24/nasa/index.html]. There are two extended “geodetic footprints” in our northern 
Kodiak geodetic network (Table 1). The KODl footprint includes a station, KODV (10 m from 


KODl), 



1990 [Ma et aL, 1990]. The footprint also includes the permanent station KODK (located within 1 km 
of KODl) constructed in late 1998 as part of an educational and outreach program with NASA’s 
Goddard Space Flight Center, Kodiak Island High School, and Thales Navigation. The second extended 
footprint includes our station MILB located near the U.S. Coast Guard station KODl, installed in 1997 
as a Continuously Operating Reference Station (CORS) of the National Geodetic Survey. Additional 
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nearby stations within the two footprints were surveyed once or twice for back-up in case the primary 
site was destroyed (Table 1). 

The U.S. Geological Survey (USGS) surveyed three stations (KRLK, AHKI, SITK) in the 
southern part of Kodiak Island as well as additional stations on the Alaska Peninsula in 1993, 1995, and 
1997, occupying each station for at least 48 hours [, Savage et aL, 1999]. The southern Kodiak stations 
were surveyed again in 2001 by the USGS and the University of Alaska. However, since there were 
large nearby earthquakes between 1999 and 2001, we used only the Kodiak Island observations through 
1997 in our analysis. 

We obtained the velocities of the GPS stations using the GAMIT/GLOBK software [King and 
Bock, 2003; Herring, 2003] and the procedure described by McClusky et aL, [2000], combining the 
observations from our field receivers with those of the global IGS network processed at the Scripps 
Orbit and Permanent Array Center (SOPAC) [Bock et aL, 1997]. To define a North America reference 
frame, we minimized the horizontal velocities of 11 stations within the non-deforming regions of the 
North America and Pacific plates (Table 2). The residual motion of these 1 1 stations is 0.8 mm/yr. 

In order to evaluate the uncertainties of our velocity estimates, we computed and examined time 
series of position within each survey and over the full span of our data for stations with at least 3 
observation time periods (Figure 3). With the a priori uncertainties we assigned to the phase 


-observations,- the ; 


we have found from analysis of similar but more extensive data sets to result in realistic uncertainties in 


velocities estimated from two or more years of observations (see, e.g., McClusky et aL, 2000). For those 
stations in the Kodiak network whose unsealed nrms scatters were greater than 1.0 we down-weighted 


the position estimates to achieve values closer to 0.7. Our estimates of vertical motions are less reliable 
than those for horizontal motion in part because of larger uncertainties in the estimates of height from 
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the short observation scenarios used for many of the surveys and in part because of larger vertical 
uncertainties in reference frame. 

Velocity results 

The estimated velocities in a North America fixed reference frame are given in Table 3 and the 
well-determined horizontal velocities are shown in Figure 2. The station velocities are plotted as treneh- 
nonnal and trench-parallel components in Figure 4 and listed in Table 3. For both the northern and 
southern Kodiak stations, the trench-normal velocities decrease nearly uniformly with distance from the 
trench (Figure 4)., 

We have a small number of stations located within 1 km of the extended footprints to assess the , 
internal consistency of the velocities.. The stations within the individual KOD1 and KODI footprints 
have similar horizontal velocities within the la level of uncertainty (Table 1); the vertical uplift rates, 
however, show greater variability. From a solution in which the estimated vertical rates of the nearest 
North American stations (at Penticton, B.C., and Fairbanks) are less than 3 mm/yr, we estimated uplift 
rates of 12-18 mm/yr for stations between 170 and 230 km from the trench. There is 10 mm/yr of 
differential uplift between CL AM/KODI/KOD K/S KIO/KRLK and KOD l/MILB/SITK The individual 
station vertical results are included for completeness. Due to the large variability in these results, as a 
function of distance from the trench, we are more confident of the horizontal velocities. 



The map view of the horizontal velocities given in Figure 2 illustrates the counterclockwise 
rotation of the velocity vectors for stations in the west relative to the eastern coastal stations. If we 
hypothesize that the orientation of the down going slab is primarily controlled by the orientation of 
relative motion between the Pacific and North American plates, we would expect the velocities above 
the locked portion of the plate interface to be at approximately the N22°W orientation. To test this, we 


7 


rotated the velocities of Table 1 into a plate motion parallel component and a plate perpendicular 
component (Table 3). Most of the Kodiak stations indicate a plate parallel component of 4-6 mm/ yr 

An alternate factor that may control the orientation of the down going slab is the geometry of the 
Aleutian trench. The orientation of down going Pacific plate motion would lead to compressional strain 
at N30°W; further i n land the trench-parallel strain would be relieved as left-lateral slip on upper plate 
faults- We resolved the data into trench perpendicular (N30°W) and trench parallel components 
(N60°W) (Table 2 and Figure 4). The Kodiak Gulf of Alaska coastal stations (SITK, PAS A, KOD1) 
have orientations that are close to the trench-normal (N30°W) but MTT.B indicates a velocity at an 
orientation (N20.8°W ± 5-5°) that is closer to that predicted for the relative plate motion (~N22°W). The 
trench parallel component increases with distance from the eastern coast. These results suggest 
distributed shear strain accumulation that could be relieved on individual faults or distributed over a 
broad region. The Kodiak Island and Narrow Cape faults are the largest mapped trench parallel faults in 
this region. Mapping and paleoseismic studies of the Narrow Cape fault show it has produced 
predominately left- lateral strike-slip displacement during the Holocene [Carver et al, 2003]. 

Modeling of Kodiak Island Geodetic Observations 

Although the epicenter for the 1964 earthquake was in the Prince Williams Sound area (between 
Anchorage and Valdez) of south central Alaska, the aftershock area extended -300 km to the east and 
-800 km to the southwest (to south of the Kodiak Islands). Two areas of high slip in the 1964 event 
correspond to seismologically determined areas of high moment release: (1) the Prince William Sound 
asperity had an average slip of 18 m and the Kodiak asperity had an average slip of 10 m [Johnson et al., 
1996; Christensen et al, 1994; Holdahl and Sauber, 1994], In our discussion of the factors that 
contribute to ongoing crustal deformation across Kodiak we consider only the contribution of the 
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southern “Kodiak asperity”. Wells et at, [2003] noted that areas of high coseismic slip commonly 
occurred beneath the prominent gravity lows outlining a deep-sea terrace low. The location of the 
Kodiak asperity does not seem to be associated with a prominent forearc basin gravity low. 

In addition to the trench parallel variability in the coseismic slip in the 1964 earthquake noted 
above, changes in the coseismic slip and interseismic locking depth as a function of distance from the 
Aleutian Trench occur as well In general, the seismogenic zone along the subduction plate interface is 
bounded by the upper and lower stability transition depths [e.g. Scholz, 1990]. The seaward updip limit 
is important for tsunami generation and the downdip limit is important for seismic hazard because it is 
used to infer the landward limit of the seismic source zone (e.g. Oleskevich et al, 1999). Since the 
boundary between the region of coseismic uplift and subsidence in the 1964 earthquake occurred on the 
eastern portion of the island, earlier studies suggested that Kodiak Island is above the transition zone 
between interface slip in large earthquakes, post-seismic creep, and aseismic creep at depth [Ma et al., 
1990; Gilpin et aL, 1994; Gilpin, 1995; Stockman et al., 1997; Savage et al., 1999; Zweck et al., 2002). 
These earlier Kodiak geodetic observations and modeling, however, hypothesize different downdip 
limits for coseismic slip and interseismic locking somewhere beneath Kodiak Island. 

We used aftershocks and tsunami studies of the 1964 earthquake, as well as thermomechanical 
models to give us an approximation of the updip limit of significant coseismic slip. Aftershocks 

from the trench (toward northern Kodiak) [Algermessen et al, 1969]. Oleskevich et al., [1999] 
estimated the 100-150°C temperature isotherm limit for clay dehydration between the Kenai Peninsula 
and Kodiak to be 80-160 km landward from the trench. In the northern Kodiak region there are some 
large earthquakes near the 1600 m contour but large earthquakes have not occurred near the trench over 
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the last 20+ years (Figure 2). In contrast, in southern Kodiak, moderate to large earthquakes have 
occurred between the trench and the Alaska Peninsula (Figure 2). 

The depth of the Aleutian trench ranges from about 4000 m near the northeastern part of the 
Kodiak Islands to 5200 m a little farther south [Beikman, 1980]. The depth and curvature of the dipping 
slab in our finite element model matches the dip of the upper slab surface inferred from background 
seismicity in the region of Kodiak Island. However, the dip in our finite element model, near the trench, 
may be too steep. The dip of the shallow portion of the slab used in the finite element model is about as 
shallow as possible given the necessity to avoid an extreme aspect ratio of the elements that would 
create a numerical instability. 

< In addition to the megathrust fault (down-going Pacific plate), Holocene surface faulting in and 
near the northeastern portion of the Kodiak Islands has been mapped offshore on the Kodiak Shelf fault 
system [Plafker et al, 1994] as well as within the Narrow Cape region [Carver et al, 2000, 2003] near 
the station PASA (Figure 1 and 2). The Border Ranges fault traverses on-shore along the northwest 
coastline of the Kodiak Islands (Figures 1 and 2) but no clear evidence of Holocene slip on this fault 
system has been observed [Gilpin, 1995]. We assumed motion along the plate interface is the 

primary source of interseismic strain but we considered slip on the upper plate faults in our 
interpretation of horizontal velocities. 


— Gver-the time- and spatial -scale of- the -surfaee- -deformation measurements^- tbis^ study, the 
response of the Earth to tectonic loading depends on the rheological structure of the crust and upper 
mantle. The time-dependent response of the Earth has been shown to be important in Alaska for 
modeling geodetic data [Brown et al., 1977; Wahr and Weiss, 1980; Sauber et al, 1993; Savage and 
Plafker, 1991; Cohen, 1996; Taylor et al, 1996; Zheng et al, 1996; Freymueller et al., 2000; Sauber et 
al, 2000, Zweck et al., 2000; Freymueller et al., 2001]. Two time-dependent processes are particularly 
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important in this region of large and great earthquakes. These are fault creep occurring down-dip of the 
coseismic rapture plane and viscoelastic shear flow in the ductile portions of the lower crust and upper 
mantle [see summary in Cohen, 1999]. Both of these mechanisms are stimulated by the coseismic 
transfer of stress from the shallow seismogenic portions of the Earth to greater depths. To evaluate 
alternate mterseismic and post-seismic mechanisms that could account for our geodetic observations 
between 1993 and 2001, we calculated the predicted average surface displacement rate for alternate 
‘ Earth models. We used primarily a two-dimensional plain strain finite element grid in TECTON 
[ Melosh and Raefsky, 1981] for these calculations. The finite element grid across the plate boundary 
includes a shallow dipping subducting slab and both an oceanic crust-mantle and a continental crust- 
mantle (Figure 5a). Because of our uncertain knowledge of the Earth’s response to tectonic loading on 

the time scale of years, we explored alternate interface slip distributions and rheological models for die 

1 

post- 1964 coseismic and mterseismic time periods. Due to the low viscosities in the lower crust and 
upper mantle (~10 19 ' 20 Pa s) inferred from most studies in Alaska, we assumed the current response to 
major late Pleistocene deglaciation and sea level rise on Kodiak Island to be insignificant Recent 
glacial fluctuations on Kodiak Island are s mall and are not modeled as well. 

Simplifying the formulation given by Cohen [1999], we calculated the surface velocity over our 
time frame of interest due to the viscoelastic response to three major conceptual elements: (1) plate 
mo tio n over -fee last -30+ years, ^2) the 1964 earthquake and (3)down-dip cr^pbefow the region OT “ 
maximum slip in the 1964 earthquake. The strain accumulation associated with the first term (1) was 
calculated as the viscoelastic response to back-slip on the shallow portion of the plate interface. The 
second term is the viscoelastic response to coseismic slip in the 1964 earthquake. Although the 
coseismic term is dependent on the (coseismic) slip distribution assumed as well as the viscosity 
structure, in general this term is smaller than (1), especially 30+ years after the earthquake. The third 


11 



term, creep on the plate interface, can be time-dependent as well Studies from the Prince William 
Sound area have suggested this is an important mechanism [Brown et al, 1977; Cohen, 1996; Cohen 
Cohen et al., 1995; Freymueller et al., 2000; Cohen and Freymueller, 2001]. Although we explore a 
range of creep models, greater spatial and temporal resolution of surface velocities is needed to resolve 
the time-dependent changes in down-dip creep. 

Coseismic slip in the 1964 earthquake 

Since the occurrence of the 1964 earthquake, numerous models of slip in the 1964 earthquake 
have been proposed and we will not attempt to summarize all the studies here. As mentioned earlier, the 
more recent models include a region of high slip referred to as the “Kodiak asperity” but studies differ in 
the location and magnitude of slip in this asperity [ Christensen and Beck, 1994; Holdahl and Sauber, 
1994; Gilpin et al., 1994; Gilpin, 1995; Johnson et al., 1996]. Gilpin [1995] included the most extensive 
Kodiak tide gauge and geologic data set and he inverted for slip across three different regions of Kodiak 
Island but not other regions (Le. the Prince William Sound region that slipped in the 1964 earthquake). 
Gilpin’s results also suggest that slip offshore extending to below northern Kodiak Island was larger 
than in southern Kodiak. 

The geologic and geodetic measurements used as constraints in the recent coseismic models 
actually bracket up to 14 months of slip. Thus estimates for “coseismic slip” include short-term post- 
seismic processes such as fault slip in aftershocks, both on the plate interlace anffwithin thesurro unding 
region, down-dip creep, and any very short-term anelastic response of the Earth. Also, it is important to 
note that for most of the coseismic slip models the vertical displacements (from offset of barnacles, tide 
gauge data, leveling) were the primary constraint. For more recent data used to estimate interseismic 
strain rates, such as our geodetic observations, the horizontal rate of deformation is the better 
constrained component. 
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In this section we refine the coseismic slip models suggested by earlier studies to make our two- 
dimensional (2-D) models specific to the two regions with geodetic observations. We used the tide 
gauge data summarized in Gilpin [1995] for the “northern” and “southern” Kodiak. The 1964 coseismic 
subsidence and uplift values for the north (black, in Figure 6) and south (red) are clearly different. For 
Model 1 (MD1) we used a coseismic slip of 10 m for the depth range between 11 and 22 km (modified 
from Holdahl and Sauber, 1994; Johnson et al , 1996). We assumed a linear transition to zero slip 
between 9 and 1 1 km and between 22 and 30 km (Figures 5a and 5b, solid line without creep between 
30 and 90 km). In such a model, slip in the 1964 earthquake was shallow and any post- seismic response 
was associated with a fairly shallow (< 40 km) low viscosity zone. The second model MD2 is partially 
based on the maximum slip in northern Kodiak given in Gilpin [1995] (Figures 5a and 5b); the- 
coseismic slip is higher (MD1 above multiplied by 1.4 m) and we assumed short-term post-seismic 
creep below the seismogenic zone (up to 3 m between depths of 22 and 90 kms). In this model we 
assumed a deeper (> 90 km) low viscosity zone. For southern Kodiak we multiplied both MD1 and 
MD2 models by 0.5 (half the coseismic slip as suggested by earlier studies). The predicted coseismic 
displacements due to these alternate coseismic slip distributions are compared in Figure 6. Although 
the two coseismic models would place different boundary conditions on the top of the asthenosphere 
right after an earthquake, the available vertical displacements at the surface can not distinguishthem. 

Interseismic model (1993-2001 ) 

As discussed above, a number of mechanisms could influence the surface deformation rates 
obtained in this study, -30-40 years following the 1964 earthquake. There are numerous studies of the 
post- seismic response to the 1964 earthquake especially in the Kenai region ( Cohen et al., 1995 and see 
summary in Freymueller et aL, 2000). Additionally, in the Kodiak region, Gilpin et al., [1994] and 
Gilpin [1995] used permanent and temporary tide gauge data from the islands to estimate a down-dip. 
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post-seismic creep of approximately 3 meters during the 25 years following the earthq uake This is 
about the amount of creep inferred for the Kenai region over a s imilar time period [ Cohen et al, 1995; 
Freymueller et al, 2000]. Based on their analysis of tide gauge data. Savage and Plafker [1991] 
postulated t h a t in addition to down-dip creep, a slower viscoelastic relaxation of the low viscosity 
mantle could be important over a longer time scale. We included the viscoelastic response to 1964 slip 
in both end member models, albeit at different depths. 

Elastic dislocation model 

Locking of the plate interface at shallow depths is assumed to exert the greatest influence on the 
interseismic rate of deformation oyer 30 years after the 1964 earthquake. For a point of reference we 
present the elastic dislocation model used by Savage et al., [1999] to represent back- slip at the relative 
plate rate on a shallow main thrust zone [Freund and Barnett, 1976; Savage, 1983]. The elastic slip 
model is shown Figure 5c and the comparison to data is shown in Figure 7a,b. Savage et al [1999] used 
a nonlinear inversion program to determine the parameters for a dislocation model which best fits their 
southern Kodiak data They assumed a dip of 5° along the locked plate interface. The locked zone 
width of approximately 152 km (57 km to 209 km downdip. Figure 3 Savage et al., 1999) and a back- 
slip rate of about 57 mm/yr were estimated. Their locked zone was between depths of 5 and 18 km 

r 

(Figure 5c). If the plastic profile is shifted relative to the trench (t he assumed zero depthis further east 
due to, for instance, a shallower dip) the elastic model will fit the data better near the locked main thrust 
zone. The elastic model however does not account for data further inland. The predicted trench 
normal rate of deformation from northern and southern Kodiak shown in Figure 7a illustrates the 
important contribution the shallow locked zone makes to ongoing strain accumulation but the 
data/model discrepancy also suggests that other processes could contribute. 


Predicted time.- dependent surface displacement rate 

For the two MD1 models, we consider two processes, the elastic response to ongoing plate 
motion and the viscoelastic response to coseismic slip in the 1964 earthquake (see summary in Table 4). 
In this model the lower crust and/or upper mantle is assumed to respond viscoelastically and most of the 
elastic layer ruptured during the 1964 earthquake. In MDla we used the coseismic region as a proxy 
for the region that is locked during the interseismic time period. Specifically, between 1 1 and 22 km we 
imposed 57 mm/yr of back-slip and between 9 and 1 1 km and 22 and 30 km, the imposed plate rate of 
back-slip was reduced to zero. For model MDlb, we used a shorter locked zone that is closer in extent 
to the elastic model of Savage. In the MD1 models, the transient component of surface deformation is 
primarily due to the time-dependent viscous relaxation of a shallow (> 38 km) low viscosity zone 
instead of time-dependent, down-dip creep. For our reference model given in the figures we assumed a 
viscosity of q of 5 x 10 19 Pa s and rigidity (G) of 0.6 x 10 19 Pa. This corresponds to a Maxwell time, r, 
of - 26 years, where t = q/G. The influence of the assumed viscosity q is discussed later. Post-seismic 
relaxation of the 1964 displacement would cause subsidence and horizontal motion away from the 
trench above the 1964 rupture. The predicted horizontal and vertical displacement rates for MDla,b are 
compared to the horizontal and vertical data and the elastic model in Figures 7a and 7b. The MDla,b 
and the elastic model predictions are similar to the observed horizontal data near the locked main thrust 
zone but trench directed motion on the Alaska Peninsula is not predicted (1993-2001, 
http://quake.wr.usgs.gov/reseaich/deformation/gps/auto/KatmaiKodiak) . Future measurements on the 
Alaska Peninsula and along the western coast of Kodiak would provide more robust constraints. 

In the MD2 suite of models the low viscosity layer is deeper (> 90 km) and post-seismic creep 
from the base of the seismogenic zone down to 90 kms is included and evaluated (Figure 5a, d and Table 
4). In MD2a, b (Table 4) we assumed that the region of high coseismic slips was also the region that 
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was locked during the interseismic time period (Le. the MD2 coseismic model). For the MD2c, d we 
assumed the region of full interseismic coupling was slightly smaller than the region of high coseismic 


slip in the 1964 earthquake; this is, closer to the locked region of Savage et al, [1999] but the MD2c, d 
viscoelastic model includes down-dip creep as well We calculated the response to 10-60 mm/yr of 
creep between the down-dip extent of the locked zone to 90 km A model with little (MD2a) or no 
creep between 30 and 90 km predicts a higher horizontal rate of deformation (arc directed) than 
observed in northern Kodiak (Figure 8a). At the other extreme, a model (not shown) that includes 60 
mm/yr of creep under predicts the horizontal rate of deformation directed arc- ward. Since the amount of 
coseismic slip is different between the north and south we included the region-specific postseismic 
relaxation due to coseismic slip for the calculations; thus the post-seismic viscous relaxation is 


hypothesized to be smaller in southern Kodiak. Additionally we hypothesize that the magnitude of 
creep scales as a function of coseismic slip as well MD2c in red represents the predicted horizontal 
rates for southern Kodiak (Figure 8a). For the model that best predicts the northern horizontal data 
(MD2c) we plotted the individual contribution from each of the deformation mechanisms (Figure 8c). 
The vertical data are best predicted by the MD2 a,b models with the wider locked main thrust zone. 

Each of these models predicts different horizontal velocities for the Alaska Peninsula region 
(Figure 8a) as well The stations from the Alaska Peninsula in the Katmai volcano region have small, 

ra North: America fixed reference 
frame (1993-2001, http://quake.wr.usgs.gov/research/deformation/gps/auto/KatmaiKodiak). These 
results are similar to that observed in the Kenai which have been attributed to dip-slip creep on a portion 
of the plate interface below the locked main thrust zone [Freymueller et al., 2000]. The higher creep 
models (MD2b,d) however, over predict the amount of trench directed motion on the Alaska Pe ninsula 


Variable viscosity values and the preferred model 
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The post-seismic component is small relative to interseismic strain and creep but we explored 
various viscosity values to test the uniqueness of the viscosity used to calculate the postseismic response 
to the 1964 earthq uake . We varied the lower crust/upper mantle viscosity between 10 19 and 5 x 10 20 Pa 
s. With a viscosity of 10 19 Pa s the Maxwell time, t, is just 5 years and for 5 x lO 20 Pa s, x = 264 years. 
The models with the viscosity values of 10 19 and 5 x 10 20 Pa s do not approximate the northern Kodiak 
data. The 10 20 Pa s model approximates the general drop off in velocity but under predicts the rate as a 
function of distance from the trench. 

Discussion 

Sources of time-dependent deformation 

Based on our modeling of the geodetic observations, time-dependent creep and possibly post- 
seismic viscoelastic relaxation may be a significant source of the measured surface velocities. Due to 
the large uncertainties in our short-term vertical results, we evaluated this conjecture by e x a m i n i n g other 
data from the region of the 1964 earthquake. Long-term and temporary tide gauge measurements have 
been used as well to estimate uplift rates since the 1964 earthquake [Savage and Plafker, 1991; Gilpin, 
1995; Cohen and Freymueller, 2001]. On Kodiak temporary tide gauge measurements have been used 
to constrain coseismic slip in the 1964 earthquake and additionally these measurements were made 
between 1965 and 1993, primarily to the, west ofrour northern K odia k GPS s t atio ns [Gilpin etah, 1 994-; 
Gilpin, 1965]. These vertical measurements show less scatter, as a function of distance of from the 
trench, than our short-term vertical rates. In Figure 9 we compared the tide gauge uplift values (1965- 
1993) to the predicted uplift assuming the model parameters given by the elastic and MDla,b (a) and 
MD2b,c (b) models. The rate of uplift for the elastic model is much higher than the observed uplift 
rates for 1965-1993 and the MDla,b models predict a lower rate of uplift than is observed (Figure 9a). 
For this earlier time period, the total uplift is better predicted by the model (MD2b) with the higher 


creep rate of 25 mm/yr (Figure 9b). Other studies have suggested a time variable uplift rate due to down- 
dip creep as well Based on the tide gauge data from Woman’s Bay (< 1 km from KODK, KODI and 
the old VLBI mark), Cohen and Freymueller [2001] estimated an average uplift rate of 14 ± 1 mm/yr 
between 1967 and 1998, but they suggest there is a significant decrease in the uplift rate as a function of 
time. They estimated an uplift rate of 15 mm/yr for 1985 and 8 mm/yr for 1995. 

Recent earthquakes (southern Kodiak 1999-2001, see discussion in introduction) and possibly 
episodic creep have occurred below the region of high coseismic slip in the 1964 earthquake in the last 
decade. Around the same time as the earthquakes in southern Kodiak, the “great Alaska earthquake of 
1998-2001” occurred in the Anchorage region [Freymueller et al, 2001]. Freymueller et al. [2001] 
hypothesized that the southeast directed surface displacement was due to creep on a segment of the plate 
interface down-dip from the locked main thrust zone. 

Long-term vertical uplift rates compared to our geodetic uplift rates 

Landsat-7 and a Shuttle Radar Topography Mission (SRTM) derived digital elevation model 
(DEM) were used by Carver et al [2003] to identify a coastal marine terrace across the northern portion 
of Kodiak. The terrace includes broad planar surfaces up to several kins wide that border the modem 
coast line. The terrace increases in elevation uniformly from about 15 m on the west of the Island 
(Region IH of Figure 1) to about 40 m in the central region (Region II of Figure 1). Based on the 
assumption that the surface had been cut by shoreline planation during the last interglacial (probably an 
oxygen isotope stage 5e, 120,000 - 130,000 ka). Carver et al. [2003] estimated an average uplift rate of 
0.15 mm/yr for Region HI to 0.3 mm/yr for Region II (Figure 1). The uplift estimated from the terrace 
data suggest that within the transition region of the 1964 coseismic uplift and subsidence (Region I of 
Figure 1) the highest rate of uplift is associated with crustal shortening on upper crustal faults. For 
instance, at Narrow Cape on Kodiak’s northeast coast, the elevation of the terrace is about 80 m and it is 
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100 m on nearby Ugak Island [Carver et al, 2003]. Mapping and trench studies of several lineaments 
identified from Landsat-7 and SRTM data show evidence for late Pleistocene and Holocene surface fault 
rupture in the vicinity of the Kodiak launch facility at Narrow Cape [Carver et aL, 2003]. 

Summary 

The new geodetic results reported in this study from northern Kodiak Islands document the 
c h a n ge in deformation rates across a transition region between uplift and subsidence in the 1964 
earthquake. Earlier VLBI results indicate that a simple elastic dislocation model used to represent 
interseismic strain accumulation was inadequate to account for the observed horizontal and vertical rates 
as a function of distance from the trench [Ma et aL, 1990]. Later studies further indicate the importance 
of time-dependent uplift rates following the 1964 earthquake [Savage and Plafker, 1991; Cohen and 
Freymueller, 2001]. Permanent GPS stations across Kodiak Island and on the Alaska Peninsula could 
provide important constraints on down-dip, time-dependent processes. In this study, we suggest that the 
recent horizontal (trench perpendicular component) and earlier vertical uplift data across Kodiak Island 
can be accounted for by the viscoelastic response to back-slip representing a locked main thrust zone, 
down-dip creep, and coseismic slip in the 1964 earthquake. Northern and southern Kodiak have 
different interface slip histories and interseismic strain predictions. A change in the orientation of the 
ho rizontal vel ocit y vectors occurs above, the down-hip, segment of the locked maim thrust -zone. We 
hypothesize that the horizontal, trench-parallel component of slip will be released as left-lateral slip on 
trench parallel, upper plate faults. 
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Figure Captions 

Figure 1. Geologic map of Kodiak Islands region of Alaska (modified from H_ M. Beikman 
1980). The location of the northern Kodiak stations are given by circles and the southern Kodiak 
Island stations in the Katrnai network of Savage et al. [1999] are given by squares. The location 
of Kodiak Island (KI) on an outline of the state of Alaska is given in the upper right hand comer. 
The scale on the lower right is given in kilometers. Region L Eastern coastal regions includes 
Narrow Cape fold and thrust belt, the Narrow Cape fault, and along the western boundary of 
Region I, the Kodiak Island fault. Brown (IT) = Lower Tertiary Ghost Rock Formation, It 
Tan(Tm) = Miocene Rocks, Narrow Cape formation. Yellow (To) = Oligocene stratified 
sedimentary sequence; Region II: Late Cretaceous Belt. This central region appears to have 
uplifted more rapidly than the region northwest of the Border Ranges fault resulting in a 
northwest tilt of the Islands [Gilpin , 1995], bright Green (uK) = upper Cretaceous stratified 
sedimentary sequence, Kodiak formation. Pink (Txif) = Paleocene felsic, intrusive rocks, 
Kodiak Batholith; Region HI: This region of volcanic arc basement forms a backstop northwest 
of the Border Ranges fault (BR), Green (Kj3) = Lower Cretaceous and Upper Jurassic rocks, 
Uyak formation, Maroon(u) =Precambrian ultramafic, Turquoise,(uTr) = Triassic volcanics, 
Shuyak formation. Pink with white (Jii)= Jurassic Blueschist-Greenschist, Raspberry Schist. BR 



northwest across the Island [Gilpin, 1995]. 

Figure 2. Map showing representative GPS station velocities from the northern Kodiak network 
(Table 1) and the southern Kodiak Island stations from the Katrnai network [Savage et al, 1999] 
with our computed horizontal velocity and 95% confidence error ellipses (in red). Focal 
mechanisms derived from Harvard Centroid Moment Tensor and USGS Fast Moment Tensor 
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solutions (1977-2004, Dziewonski and Woodhouse 1983; Ekstrom, 1994; 
httpj7neic.usgs.gov/neis/sopar) were plotted using the program iGMT [Becker and Braun, 1998], 
The event location, fault plane parameters, and date are given in Table 5. Note: The longitude 
is given in °E. 

Figure 3. The detrended north, east, and vertical time series of station position for stations from 
the Kodiak GPS northern and southern networks that have been observed at least 3 times. 

Figure 4. The station velocities resolved into trench-normal and trench-parallel components as a 
function of distance from the trench (taken to be from the southeastern most extent of the 4000 
or 5800 m bathymetry contour). The error bars are two-sigma. The station names are indicated 
along the top of the graph, the black data points are from northern Kodiak and southern Kodiak 
stations are given in red. Positive velocities for the trench normal component indicate away 
from the trench. Negative velocities for the trench parallel component indicate a west directed 
velocity. 

Figure 5 (a) Two dimensional finite element grid used in MD1 and MD2 models (modified from 
Cohen [ 1996], with the details of the model geometry given in Table 3 in that reference). The 
“zero slab depth” location in our finite element and dislocation models is seaward (primarily 
east) of the trench and corresponds to the projection of the slab to the surface. A low viscosity 
zcneuccu rs belo w -a-depth of 3&dma-ifl J MBP andbelow90 kminthe MD2 models.' In our - 
reference model a viscosity of 5 x 10 19 Pa s is assumed; this corresponds to a Maxwell time, x, 
equal to 26 years; this is short relative to the estimated repeat time of major earthquakes in the 
northern Kodiak region of 300-500 years [Gilpin, 1995]. The dark line between the downgoing 
Pacific plate (in red) and the overriding continental plate (in green) represents the region of 
reduced slip during the 1964 earthquake and partial locking during the interseismic time interval 
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in some models. Within the "Key" are the assumed elastic parameters (E) and viscosity p for the 
reference model. 

b. Coseismic slip (meters) in the 1964 earthquake as a function of depth for the MD1 (solid line) 
and MD2 (dashed line) models. The MD2 model includes reduced slip between 22 and 90. 

c. Input slip (mm/yr) for interseismic time period as a function of depth for Elastic (dotted), 
MDla (dashed) and MDlb (solid line). 

d. Input slip (mm/yr) for interseismic time period as a function of depth for MD2a (dashed), 
MD2b (dotted), and MD2c (solid) and MD2d (dash-dot). 

Figure 6. Predicted coseismic vertical displacement as a function of distance for the two slip 
models in Figure 5b along with coseismic data (asterick with error bar for the north and diamond 
with error bar in red for the south) from Plafker [1969). Model 1 (MD1) is given by a solid line 
and Model 2 (MD2) by a dotted line, in black for northern Kodiak and in red for southern 
Kodiak. The coseismic data of Plafker [1969] given here corresponds to a time interval of up to 
14 months (see text). The down-dip slip in MD2 between 22 and 90 km is assumed to represent 
a longer post-seismic time interval 

Figure 7. (a) Comparison of the geodetieally estimated rate of deformation (trench normal) to the 
predicted horizontal velocity from the Elastic (dotted. Savage et al, [1999]) and the MDla 
- fsehdF MDla(sehd-red} for- southern Kodiak and MB lb (dashed) viscoelastic models. The 
northern Kodiak trench-normal data is given in black and the southern Kodiak data in red. Note 
that none of these models predict trench directed velocities at distances of -300-400 km (Alaska 
Peninsula), (b) same models as given in a. but for the vertical displacement rate as a function of 
distance from the trench. The vertical uplift rate at KRLK is not matched by any of the models. 
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Figure 8. (a) Comparison of the predicted horizontal velocities for the four MD2 models of 
Table 4, plus southern Kodiak for the preferred model of MD2c (red), to the observed horizontal 
velocities, (b) the observed versus vertical displacement rate for the four MD2 models, plus 
southern Kodiak model of MD2c (red), (c) the observed versus predicted horizontal velocities 
for the MD2c model including the individual elements (dotted=interseismic, 
dashed=postseismic, and dash-dot=creep) and all three components together (total solid Iine). 
Data description as given in Figure 7. And (d) the predicted horizontal velocity for the MD2c 
model assuming four different viscosities 10 19 , 5 x 10 19 , 10 20 and 5 x 10 20 Pa s and for southern 
Kodiak MD2c (solid red) for a viscosity of 5 x 10 19 Pa s. 

Figure 9. Comparison of the predicted vertical displacement between 1965 and 1993 assuming 
the model parameters for MD2b (dashed line), MD2c (solid line), MD2c (solid red) over a 28 
year time period compared to the temporary tide gauge data of Gilpin [1995]. The higher creep 
rate model of MD2b better matches the uplift rate observed 1965-1993. 
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Table 2. North America fixed reference frame stations. 


Station North Rate 
(mm/yr) 

N 0 

(mm/yr) 

East Rate 
(mm/yr) 

E 0 

(mm/yr) 

Vertical 

(mm/yr) 

v a 

(mm/yr) 

North American stations 





STJO 

0.82 

0.86 

-0.22 

0.87 

-0.28 

0.87 

BRMU 

0.70 

0.62 

0.15 

0.63 

-1.42 

0.85 

ALGO 

-0.33 

0.63 

0.25 

0.64 

2.98 

0.73 

NUB 

0.50 

0.32 

0.53 

0.34 

-1.14 

0.82 

MDOl 

0.11 

0.30 

-0.37 

0.33 

0.23 

0.82 

PIE1 

-0.25 

0.78 

-0.26 

0.80 

2.08 

0.80 

YELL 

-0.44 

0.25 

0.38 

0.25 

3.86 

1.13 

WILL 

0.78 

0.98 

-0.08 

1.00 

-0.08 

1.27 


Pacific Plate Stations 


CHAT52.87 

0.89 

-38.22 

0.88 

1.14 

1.30 


KWJ1 47.17 

1.13 

-69.44 

1.33 

-3.76 

2.86 


MKEA 51.62 

0.87 

-59.75 

0.91 

-2.67 

2.83 



Note: The correlation between components is less than 0.03. 
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Table 5 Event locations, fault plane parameters, and date. 


Event 

Num. 

Lon 

(°E) 

Lat 

(°N) 

Strl 

Dipl 

Rakel Str2 

Dip2 

Rake 2 

1 

-154.43 

58.09 

42 

65 

22 

302 

70 

153 

2 

-151.21 

56.30 

265 

8 

119 

55 

83 

86 

3 

-154.16 

56.20 

244 

13 

94 

60 

77 

89 

4 

-151.38 

56.91 

215 

9 

78 

47 

82 

92 

5 

-155.43 

56.95 

85 

59 

-18 

185 

75 

-147 

6 

-154.18 

58.62 

334 

25 

-164 

229 

83 

-65 

7 

-154.06 

58.58 

319 

50 

168 

56 

81 

41 

8 

-153.38 

56.11 

220 

18 

75 

56 

73 

95 

9 

-153.20 

58.88 

89 

64 

12 

353 

80 

153 

10 

-153.21 

56.06 

241 

11 

105 

46 

79 

87 

11 

-152.28 

56.30 

228 

8 

85 

53 

82 

91 

12 

-154.64 

57.70 

319 

28 

-176 

225 

88 

-62 

13 

-153.17 

58.30 

176 

62 

-154 

73 

67 

-31 

14 

-153.09 

56.89 

212 

9 

61 

62 

82 

94 

15 

-155.48 

56.34 

0 

70 

173 

93 

84 

20 

16 

-154.87 

57.24 

358 

29 

-148 

239 

75 

-65 

17 

-153.81 

56.80 

61 

33 

-67 

214 

60 

-105 

18 

-149.96 

58.49 

141 

23 

-150 

23 

79 

-70 

19 

-152.83 

57.03 

58 

4 

-66 

214 

87 

-91 

20 

-151.54 

56.14 

218 

20 

60 

70 

72 

101 

21 

-152.07 

56.73 

299 

5 

159 

50 

88 

85 

22 

-152.86 

56.48 

209 

11 

101 

18 

79 

88 

23 

-152.58 

56.21 

225 

9 

74 

61 

81 

93 

24 

-154.78 

58.15 

64 

49 

-17 

165 

78 

-138 

25 

-154.35 

57.35 

357 

63 

-180 

267 

90 

-27 

26 

-154.38 

57.22 

350 

37 

-161 

245 

79 

-55 

27 

-154.04 

57.24 

5 

75 

178 

96 

88 

15 

28 

-149.29 

57.47 

302 

55 

136 

61 

55 

45 

29 

-154.44 

57.41 

350 

10 

-139 

219 

84 

-83 
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15 

-96 

Q1 9 

zrxrz 7 

75- 

- 88 - 

31 

-152.84 

58.97 

330 

56 

-174 

236 

85 

-35 

32 

-150.18 

58.65 

82 

63 

24 

341 

69 

150 

33 

-154.22 

57.54 

356 

47 

-161 

253 

76 

-44 

34 

-154.24 

57.58 

224 

42 

-36 

341 

67 

-127 

35 

-153.70 

58.00 

98 

60 

35 

349 

60 

145 

36 

-153.56 

56.99 

224 

8 

74 

60 

82 

92 

37 

-153.63 

56.51 

225 

7 

92 

43 

83 

90 

38 

-150.99 

57.06 

332 

66 

170 

66 

81 

24 

39 

-154.97 

58.82 

63 

64 

26 

321 

67 

151 

40 

-155.11 

57.63 

278 

35 

73 

119 

57 

102 

41 

-151.26 

58.50 

218 

32 

-31 

335 

74 

-118 

42 

-151.44 

58.33 

154 

8 

-14 

258 

88 

-98 


Date 

(MMDDYR) 

112777 
071978 
070680 
090682 
112683 
032384 
102785 
030686 
052386 
091686 
052087 
061689 
030890 
052990 
080192 
091292 
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121594 
030396 
120896 
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030800 
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Popular Summary of “Crustal Deformation and the Seismic Cycle across the Kodiak 
Islands, Alaska” by J. Sauber(NASA GSFC), G. Carver (Humboldt University), S. Cohen (NASA 
GSFC) and R. King (MIT) 

The Kodiak Islands are located approximately 120 to 250 km from the Alaska- Aleutian Trench 
where the Pacific plate is underthrusting the North American plate at a rate of about 57 mm/yr. The 
southern extent of the 1964 Prince William Sound (M w = 9.2) earthquake rupture occurred offshore and 
beneath the eastern portion of the Kodiak Islands. In this study we report global positioning system 
(GPS) restilts (1993-2001) from northern Kodiak Island that span the transition between the 1964 uplift 
region along the eastern coast and the region of coseismic subsidence further inland. We used these 
geodetic results to look at the amount of slip near Kodiak during the great 1964 earthquake and we 
examined the relation of the coseismic slip to crustal deformation measured in the 30 years following the 
great earthquake. In addition to strain accumulating that will be released in the next large subduction 
earthquake, we suggest that 4-8 mm/yr of slip is accumulating that will be released as left-lateral strike- 
slip motion across the inland faults of Kodiak Island. Based on the pre-1964 and post-1964 earthquake 
history, as well as the pattern of interseismic earthquakes across the plate boundary zone, we 
hypothesize that in southern Kodiak some strain is released in moderate to large earthquakes between 
the occurrences of great earthquakes like the 1964 event. In northern Kodiak, however, the main thrust 
zone is locked and will eventually be relieved in another large subduction zone earthquake in the future. 


